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Abstract 

An extension of the Newell- Whitehead-Segel amplitude equation covariant 
under abritrary rotations is derived systematically by the renormalization 
group method. 
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The present day theory of pattern formation to a large degree rests on the derivation and 
exploitation of amplitude equations. In particular, for the formation of 2-dimensional pat- 
terns by spontaneous symmetry breaking at a finite wavenumber in isotropic 2-dimensional 
layers the amplitude equation due to Newell and Whitehead and to Segel [ffj is of funda- 
mental importance. It is derived, e.g. for the Benard instability, close to the onset of the 
pattern, by assuming the presence of an ideal one- dimensional roll pattern and certain scal- 
ing properties of its small amplitude and its spatial and temporal variations. The small 
parameter is the difference of the bifurcation parameter from its value at onset. 

The resulting amplitude equation, despite of its great fundamental and practical impor- 
tance, has a well known short-coming: it does not respect the full rotation invariance of the 
two-dimensional system and the fact, that in principle the rotational symmetry is broken 
by the pattern spontaneously rather than by any external agent. This short-coming is a 
direct consequence of the method of derivation, which singles out a particular direction for 
the main pattern, and even more importantly, makes use of anisotropic assumptions for the 
scaling of the spatial variation of the pattern. 

Two recent developments now permit to overcome this problem 0. First, Gunaratne et 
al recently proposed an extended form of the Newell- Whitehead-Segel (N-W-S) equation 
which they demonstrated to fully respect the rotational symmetry of the original system. 
They also provided a derivation of their equation by including, in a given order in the bi- 
furcation parameter, symmetry restoring terms, which, indeed appear in higher order of the 
N-W-S scheme. While this recombining of terms from different orders of the expansion sug- 
gests that their equation is the correct symmetric extension of the N-W-S equation, before 
one can accept this as a fact, a fully systematic (if not mathematically rigorous) derivation 

2 



of this equation must be provided. Clearly, a new method, avoiding any anisotropic scaling 
assumptions, is required. Recently, such a method has become available with the renormal- 
ization group method by Chen et al applied to the derivation of amplitude equations 
|5||J. Here we shall employ this method to give a systemtic derivation of the rotationally 
symmetric extension of the N-W-S equation. Curiously, Chen et al |J and also Kunihiro || 
already used this method to rederive the N-W-S equation, seemingly without invoking the 
anisotropic scaling assumption of the original derivation. However, as we shall see, such scal- 
ing assumptions in fact crept inadvertently into their derivation by the assignement which 
out of a number of space-time secular terms would be 'most singular'. Avoiding the scaling 
assumption implicit in their choice and simply keeping all singular terms, but otherwise 
following closely the treatment in |5j we are able to achieve our goal. 

For concreteness and simplicity, we shall present the derivation starting from the two- 
dimensional Swift-Hohenberg equation [^] 

where e is the bifurcation parameter, and if) was scaled to make the i[) and -0 3 -terms both of 
the same order in e. 

The derivation starts with the ansatz 

ip = 2Re{Ae ikx + e^ l (x,y) + ---} (2) 

with amplitude A which is constant to lowest order in e. The ansatz (|J) breaks the rotational 
symmetry in the (x, y)-plane. However, due to the rotationally symmetric form of eq. ([[]) 
the resulting amplitude equation for A will actually be covariant under rotations of k = ke x 
in eq. @ in a sense to be made concrete below. Inserting eq. @ in eq. (H) one obtains an 
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equation 



with 



^2L-\ ipi — (1 - 3\A\ 2 )Ae tkx + non-resonant terms 

\i=l ) 



La = _ 4A; 2 (£--ik) , L, = 2 (I- - ik V 



(3) 



<9x y y<9x J dy 

Aik I zA; ) — — — , L 7 
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ydx J dy 2 ' dy A 
The Li turn out all to commute. Now, the space-time secular part of a general solution of 
eq. (H) is constructed in the form 



^ = 2Re\Ae lkx + e(l - 3\A\ 2 )Ae ikx (c x t + C 2 ^- + C 37 ^- + C A ' ' 
L \ 24 zAik 
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where we have displayed all secular terms to order e and the dots denote non-secular or 
higher order terms. The Ci are integration constants and satisfy J2%Ci = 1. Incidentally, in 
H and the terms ~ t, x 2 , xy 2 , y A were considered as 'most singular' and only those terms 
were kept there. Because this assignment implies different scalings for the other terms and 
breaks the rotational symmetry we shall not use it here. This accounts for the differences 
of our results below from those in [^]|| . 

The renormalization group method [||] now proceeds by (i) introducing an arbitrary 
'regularization point' X, Y, T and to reorganize eq. (|4]) as an expansion around this point, 
rather than the origin; (ii) splitting in the amplitudes of the secular terms x n = (x n — X n ) + 
X n and similar for y n , t n where n is any of the exponents of x, y, t appearing there; (hi) 
absorbing X n , Y n , T n , by a redefinition of the amplitude A which thereby becomes a function 



A(X, Y, T). The equation we are left with is eq. (Q) with the replacements A — > A(X, Y, T), 
t -> i - T, x n -> x n - X n , y n -> y n - F n , but e* 



The symmetry which gives rise to the renormalization group is the fact that ^ is inde- 
pendent of X, Y, T. Therefore, all derivatives of ip with respect X, Y, T must vanish. To first 
order in e one obtains the conditions 



4ik 



dA 
dT 
d 4 A 
dX* 
d 3 A 
dX* 
d 2 A 
dX 2 

o 4 a 

''8X 2 dY 2 
d 3 A 



Aik 
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dXdY 2 
d A A 



dY 4 



edA(l-3\A 
eC 2 A(l - 3\A 
eC 3 A{l - 3\A 
eC A A(l - 3\A 
eC 5 A(l - 3\A 
eC 6 A(l - 3\A 
eC 7 A(l - 3\A 



2 ) =0 
2 ) =0 
2 )=0 
2 ) =0 
2 ) =0 
2 )=0 
2 ) =0 



Adding these equations by using J2% C% — 1 an d identifying X,Y,T — > x, y, t one obtains 
ip = 2Re{A(x, y, z, t)e zkx } instead of (§) and the amplitude equation, independent of the Cj, 
to order e 



4k 2 n 2 x A 



(5) 



with the symbol O x introduced in H] as 



□x = e x ■ V -V 
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The N-W-S equation is also of the form (|5|), but there n 2 is replaced by (Jj — j^-^p) 2 - It is 
obtained for special solutions of (Bf) with Ci = C3 = C5 = 0. This completes our derivation. 



It is clear that the method is systematic and may be extended to higher order in e. It can be 
seen that the particular form D x of the spatial derivative operator is a direct consequence 
of the form of all secular terms in (f|), which are in turn a consequence of the rotational 
invariance of the derivative operator in ([I]). Hence O x will automatically be generated by 
this method starting from any rotationally invariant set of equations in place of ([]]), like e.g. 
the Boussinesq equations for thermal convection. 

The covariance of eq. (|5|) under rotations has been discussed in ||. It rests on the 
remarkable property n x e lAkxX = ||, where Ak x = k($) — ke x = fc((cos# — l)e x + sm$e y ), 
and the even more remarkable property 

-iAk x x n iAk x x _ n 
t- u xC — L - 'x' J 

where e x > = cos$e x + sin^e y . It is the second property (which is independent of the first 
because D x is of second order) which makes in eq. (|^) the transformation of the amplitude 
A — > Ae lAkx ' x equivalent to the rotation ke x — > ke x / of the wave- vector. 

As many properties of patterns like e.g. defects and textures depend in a crucial manner 
on the underlying symmetry and the fact that it is spontaneously, not externally, broken, 
the covariance of eq. (|5]) under rotations must be considered an important improvement of 
the fundamental N-W-S equation. In addition we may remark that the derivation we have 
given would still apply if noise is added on the right hand side of (0) providing a rotationally 
covariant extension of the noisy amplitude equation derived in || . 

Like the N-W-S equation the amplitude equation ([|) has a potential, i.e. it can be 
written in the form ^ = — j^- with the Lyapunov functional 
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However, unlike the potential for the N-W-S equation || (where again D x is replaced by 
(J^ — ^ J^)) the potential (|6]) has minima, the attractors of eq. (||), for A = (3) -1 / 2 e tAk ^' x 
with arbitrary angle In the noisy case W({A, A*}) = N exp[—<f)/Q] is the steady-state 
distribution of the fluctuating amplitude. Here Q is a measure of the noise intensity . For 
the case of Benard convection at finite temperature it has been calculated in ||. Again 
rolls in arbitrary directions in the (x, ?/)-plane appear with equal weight in this distribution, 
in contrast to the case of the N-W-S equation where fc-vectors parallel to the x-axis are 
strongly preferred ||. 

In summary, we have given a systematic derivation of the amplitude equation for the 
formation of roll patterns in 2-dimensional systems by spontaneous symmetry breaking. The 
derivation is free from anisotropic scaling assumptions and establishes a rotationally form- 
invariant extension of the N-W-S equation as the fundamental amplitude equation for such 
systems. 
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